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Abstract. The seasonal cycle and interannual variability of two estimates of soil (or mineral ) 
dust aerosols are compared: Advanced Very High Resolution Radiometer (AVHRR) aerosol optical 
thickness (AOT) and Total Ozone Mapping Spectrometer (TOMS) aerosol index (AI). Both data 
sets, comprising more than a decade of global, daily images, are commonly used to evaluate 
aerosol transport models. The present comparison is based upon monthly averages, constructed 
from daily images of each data set for the period between 1984 and 1990, a period that excludes 
contamination from volcanic eruptions. The comparison focuses upon the Northern Hemisphere 
subtropical Atlantic Ocean, where soil dust aerosols make the largest contribution to the aerosol 
load, and are assumed to dominate the variability of each data set. While each retrieval is sensitive 
to a different aerosol radiative property — absorption for the TOMS AI versus reflectance for the 
AVHRR AOT — the seasonal cycles of dust loading implied by each retrieval are consistent, if 
seasonal variations in the height of the aerosol layer are taken into account when interpreting the 
TOMS AI. On interannual time scales, the correlation is low at most locations. It is suggested 
that the poor interannual correlation is at least partly a consequence of data availability. When 
the monthly averages are constructed using only days common to both data sets, the correlation 
is substantially increased: this consistency suggests that both TOMS and AVHRR accurately 
measure the aerosol load in any given scene. However, the two retrievals have only a few days in 
common per month so that these restricted monthly averages have a large uncertainty. Calculations 
suggest that at least 7 to 10 daily images are needed to estimate reliably the average dust load 
during any particular month, a threshold that is rarely satisfied by the AVHRR AOT due to the 
presence of clouds in the domain. By rebinning each data set onto a coarser grid, the availability of 
the AVHRR AOT is increased during any particular month, along with its interannual correlation 
with the TOMS AI. The latter easily exceeds the sampling threshold due to its greater ability to 
infer the aerosol load in the presence of clouds. Whether the TOMS AI should be regarded as 
a more reliable indicator of interannual variability depends upon the extent of contamination by 
sub-pixel clouds. 


1. Introduction 

Despite the increase in globally averaged surface tem- 
perature over the past few decades, some regions of 
the world have actually cooled [Houghton et ai , 1996; 
Hansen et ai , 1999]. Regional variations in the observed 
warming have drawn attention to aerosols, whose distri- 
butions are highly variable in space and time [Charlson 
et ai , 1992; Duce , 1995; Tegen et ai, 1996]. Aerosols per- 


turb the Earth’s climate in two ways: first, by absorbing 
and scattering solar and terrestrial radiation (the ‘direct’ 
effect), and by acting as condensation nuclei for cloud 
droplets, thus changing cloud optical properties and life- 
times (the ‘indirect’ effect). By mass, the most common 
aerosols are soil particles lofted from dry land surfaces 
by the wind [ Andreae , 1995; Tegen et ai, 1997]. Soil 
dust (also known as mineral dust) is one of the largest 
contributors to aerosol radiative forcing, and may have 
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increased in time as a consequence of human disturbance 
of soil surfaces through agriculture, deforestation, and 
overgrazing [Tegen et a/., 1996]. In order to quantify the 
natural and anthropogenic components of radiative forc- 
ing by soil dust aerosols and their effect upon climate, 
one needs continuous and global measurements of both 
their horizontal and vertical distribution. 

Original estimates of aerosol distribution and trans- 
port were derived from ground-based measurements. For 
soil dust aerosols, the observational record at certain lo- 
cations extends back several decades to the 1960’s [ Carl- 
son and Prospero , 1972; Prospero et a/., 1981; Prospero 
and Nees , 1986; Duce , 1995]. These observations al- 
low inferences about particular source regions, but they 
are not sufficient to map the global distribution of dust 
aerosols, which has large regional variability. In the past 
two decades, satellite measurements have provided de- 
tailed inferences of the horizontal distribution of dust 
aerosols over specific regions as well as over the entire 
globe [Moulin et a/., 1997b; Husar et ai, 1997; Herman 
et a/., 1997; King et al ., 1999]. Satellite estimates of the 
seasonal cycle and interannual variability of dust aerosols 
are a useful test of models that predict the aerosol dis- 
tribution (e.g [Tegen and Fung , 1994, 1995; Tegen and 
Miller , 1998]). In particular, the observed horizontal 
distribution of soil dust aerosols can be used to evalu- 
ate atmospheric general circulation model (AGCM) pa- 
rameterizations of dust uplift, transport, and deposition. 
AGCMs are essential to calculate the climate response 
to aerosol forcing [A/i7/er and Tegen , 1998], including the 
anthropogenic component, and this estimate will depend 
upon the validity of the dust parameterization. 

There are many satellite estimates of the aerosol load 
[King et al., 1999]. However, because of their global 
coverage and multi-decadal span, two retrievals are par- 
ticularly attractive to aerosol modelers in their effort 
to evaluate the distribution of soil dust calculated by 
their models: Advanced Very High Resolution Radiome- 
ter (AVHRR) Aerosol Optical Thickness (AOT) and To- 
tal Ozone Mapping Spectrometer (TOMS) Aerosol Index 
(Al). This paper compares the seasonal cycle and inter- 
annual variability of soil dust aerosols over the Northern 
Hemisphere (NH) subtropical Atlantic as inferred from 
these two retrievals. While each dataset contains ac- 
knowledged uncertainties [Husar et a/., 1997; Herman 
et al., 1997; Torres et al., 1998], consistent features on 
various time scales could be used to evaluate dust pa- 
rameterizations in the AGCM. 

Comparison of the two retrievals requires some care. 
Both TOMS and AVHRR are sensitive to the AOT 
within the column, but the TOMS Al depends addition- 
ally upon the height of the aerosol layer, which has a 
pronounced seasonal cycle [Carlson and Prospero , 1972; 
Chiapello et ai, 1995]. In section 2, we describe each 
retrieval, and the conditions under which an AOT and 
AI can be compared. In section 3, we show that each 
retrieval implies a consistent seasonal cycle of the dust 
load, if observed seasonal variations in the height of 


aerosol transport are accounted for when interpreting the 
TOMS AI. In contrast, the agreement between interan- 
nual variations in the two retrievals — measured by cor- 
relating monthly averages of each retrieval after remov- 
ing the seasonal cycle — is much smaller. In section 4, 
we show that the interannual correlation is substantially 
increased when the monthly averages are formed using 
only days for which both retrievals are available. This 
high correlation suggests that TOMS and AVHRR form 
similar estimates of the dust load within any particular 
scene. We suggest two reasons for the poor interannual 
correlation of the monthly averages, when constructed 
using all the days available in each data set. First, the 
presence of clouds may preclude a sufficient number of 
AVHRR retrievals in any particular month to allow the 
accurate estimation of the monthly average along with 
month to month variability. Second, undetected clouds 
at the sub-pixel scale may introduce uncertainty into the 
TOMS AI. Our conclusions are presented in Section 5. 


2. Datasets 

TOMS detects the presence of absorbing aerosols such 
as soil dust, smoke, and soot, based upon departures 
of backward scattered UV radiation from pure Rayleigh 
scattering. An aerosol index (AI) is constructed using 
the spectral contrast between the 340 and 380 nm chan- 
nels [Herman et ai, 1997; Torres et ai, 1998]. The AI de- 
pends upon the aerosol optical thickness (AOT), height 
of the aerosol layer, optical properties of the aerosol par- 
ticles, and instrument viewing angle. [Hsu et ai, 1999] 
use measurements from ground-based sun photometers 
to show that daily variations of the AOT and A I have 
an approximate linear relationship; weak spatial and sea- 
sonal variability of the regression coefficient relating the 
AOT and AI is attributed to variations in the height 
of the aerosol layer. TOMS detects aerosols not only 
over the water but also over land, because the UV re- 
flectivity of the Earth’s surface is very small except for 
regions covered by snow and ice [Herman et ai, 1997]. 
However, because a low altitude absorbing aerosol layer 
creates only a small perturbation upon the upwelling 
UV, the AI cannot distinguish soil dust within roughly 
a kilometer of the surface from more weakly absorbing 
aerosol types or clouds. [Herman et ai, 1997; Torres 
et ai, 1998; Hsu et ai, 1999]. 

In contrast, AVHRR measures reflectance at the 
630 nm wavelength to derive an AOT corresponding to 
all aerosol types — including weakly absorbing aerosols 
like industrial sulfates — throughout the entire thickness 
of the atmospheric column. However, AVHRR is lim- 
ited to retrieval over the ocean because of the difficulty 
presented by the nonuniforin reflectivity of the Earth’s 
land surface at 630 nm. [Husar et ai, 1997]. AVHRR 
derives an aerosol optical thickness assuming Mie scat- 
tering by spherical particles. In fact, dust aerosols are 
non-spherical, and this assumption may introduce er- 
rors in retrieval [Mishchenko et ai , 1995]. Furthermore, 
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the AOT retrieval assumes a mean particle radius of 
0.1 /mi, a value that is appropriate for sulfate aerosols 
but smaller than typical dust particle sizes [Tegen and 
Lacis , 1996; Moulin et al, 1997c]. The Mie scatterers are 
assumed to be perfect reflectors, with a single scattering 
albedo of unity; again, this is a better approximation 
for sulfate aerosols than soil dust [Hansen et a/., 1998]. 
Lastly, the ocean albedo is assumed to be independent 
of the overpass time of day. To measure the uncertainty 
introduced by these assumptions, the original AVHRR 
AOT retrieval was compared to sun photometer mea- 
surements [ Ignatov et a/., 1995]. As a result, a second 
generation retrieval of AOT was derived, which is used 
in this study. 

In the next section, we correlate dust variability as 
inferred by the TOMS AI and AVHRR AOT. This com- 
parison is complicated by two factors. First, the re- 
trievals do not detect soil dust per se , but rather aerosol 
absorption (TOMS) and reflection (AVHRR), and thus 
potentially detect additional aerosol types. We choose 
for our domain of comparison the NH subtropical At- 
lantic Ocean, because the aerosol load in this region is 
dominated by soil dust originating from North Africa 
(e.g., [Li et a/., 1996; Tegen et al ., 1997; Chiapello et al . , 
1999]). In this region, we assume that both retrievals 
predominately record variations of soil dust, excluding 
North American industrial sulfates potentially present 
in AVHRR, and carbonaceous aerosols from southern 
African biomass burning in TOMS. 

The second complication results from the slightly dif- 
ferent physical quantities represented by each retrieval: 
as noted above, the daily AI retrievals depend upon 
aerosol height, single scatter albedo, and instrument 
viewing angle, in addition to AOT. On daily and sea- 
sonal time scales, variations of the aerosol layer height 
are large enough that the AI and AOT are not directly 
comparable. To minimize the effect of synoptic and 
mesoscale variations in the height of aerosol transport 
and its effect upon the AI, we compare monthly aver- 
ages of the two data sets, rather than the original daily 
values. (The use of monthly averages also minimizes the 
importance of the instrument viewing angle, which oth- 
erwise imposes cyclic variations upon the AI with the 6 
day overpass frequency of the satellite.) Monthly aver- 
ages of the AI remain sensitive to seasonal and interan- 
nual variations in the height of the aerosol layer. During 
the NH summer, soil dust originating from Africa tra- 
verses the Atlantic. Ocean in a layer centered a few kilo- 
meters above the surface [Carlson and Prospero , 1972]. 
In contrast, wintertime transport across the Atlantic oc- 
curs adjacent to the surface [Chiapello et a/., 1995]. To 
compare the seasonal cycle of dust loading inferred by 
TOMS and AVHRR, we will attempt to account for the 
variation in layer height using the theoretical calcula- 
tions of AI by [Torres et a/., 1998]. We are particularly 
interested in whether the TOMS AI corroborates the 
amplitude of the AVHRR seasonal cycle, since certain 
dust models underestimate this cycle downwind of the 


Sahara, compared to the AVHRR retrieval (e.g. [Tegen 
et al , 1997]). 

On interannual time scales, we assume that the AI 
and AOT are directly comparable. That is, after re- 
moval of the seasonal cycle, variations in the height of 
the aerosol layer are assumed to be of secondary impor- 
tance to the AI. An identical assumption is made by [Hsu 
et al. , 1999] who interpreted interannual changes in the 
AI as a proxy for changes in the AOT. 

The AVHRR retrievals are available on a uniform 1° 
longitude by 1° latitude grid while TOMS is on a 1.25° 
by 1° grid. Both datasets are gridded at the same lat- 
itudes; in order to match the two grids with respect to 
longitude, each TOMS point is matched with the closest 
AVHRR point. Both data sets are available on a daily 
basis, but on any given day, data is missing at certain 
locations. In the case of TOMS, sea-glint at the nadir 
point prevents retrieval over ocean locations roughly ev- 
ery 6 days, which corresponds to the overpass frequency. 
For AVHRR, this frequency is around 9 days, and re- 
trievals at any location are available roughly half of the 
time, when the location is on the side of the satellite 
away from the sun. 

In practice, the presence of clouds further limits the 
availability of each retrieval. Clouds exhibit a small de- 
parture from Rayleigh scattering at the AI wavelengths 
[ Herman et al. , 1997]. Thus, a small value of the TOMS 
AI may represent a scene where the aerosol layer is 
almost entirely obscured by clouds — and cannot be 
retrieved accurately — or else a clear scene with no 
aerosols present. To distinguish these two cases, we 
choose a reflectivity threshold to be compared to the 
TOMS 380 nm reflectivity. If the reflectivity exceeds the 
threshold, say 15 percent, then we assume that the AI 
is adulterated by the presence of clouds, and discard the 
retrieval for that particular day. We also discard values 
of the AI less than zero, since this is generally the sig- 
nature of either clouds or a layer of weakly absorbing 
aerosols like sulfates. We find that the correlation be- 
tween the TOMS AI and AVHRR AOT is insensitive to 
whether a reflectivity threshold of 15 or 25 percent is cho- 
sen. The results presented in this article are based upon 
the more stringent 15 percent threshold unless other- 
wise noted. Over the NH subtropical Atlantic, where the 
dust concentration is highest (e.g., [Moulin et a/., 1997b; 
Husar et a/., 1997; Herman et a/., 1997], the presence 
of clouds and weakly absorbing aerosols cause slightly 
more than half of the daily TOMS AI retrievals to be 
discarded during any particular month; monthly aver- 
ages are constructed on average from 10 to 15 daily re- 
trievals (Fig. la). For comparison, a reflectivity thresh- 
old of 25 percent makes available roughly 5 extra days 
per month. Fig. la represents the availability averaged 
over the entire period of comparison; in the NH sub- 
tropics, availability varies from month to month with a 
standard deviation of around 4 days. 

The presence of clouds leads to a greater reduction 
in the availability of the AVHRR daily retrievals. In 
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the AVHRR AOT retrieval algorithm, clouds are distin- 
guished from aerosols by their spectral signature, using 
a test that is effective for values of the AOT less than 2 
[Husar et al. , 1997] — in fact, we found few daily values 
of AOT within our domain that exceeded 1.2, and none 
above 1.5. Over the tropical Atlantic, daily retrievals are 
typically available fewer than 10 days per month, and 
less near the Intertropical Convergence Zone (ITCZ), as 
shown in Fig. lb. The standard deviation of month to 
month variations in availability ranges from 2 to 3 days. 
We will consider the effect of data availability upon the 
correlation of the two retrievals in Section 4. 

Sub-pixel clouds that remain undetected introduce 
uncertainty into each retrieval, causing the aerosol load 
to be underestimated by the TOMS AI and overesti- 
mated bv the AVHRR AOT [Ton'es et a/., 1998]. A 
larger uncertainty is introduced into the former as a re- 
sult of its larger instrument footprint. This effect, which 
we are unable to quantify, would reduce the correlation 
between the two retrievals. 

From the nearly two decades of retrievals available 
within each data set, we examine the seven-year period 
from 1984 to 1990, which excludes any contamination 
from the volcanic eruptions of El Chichon and Pinatubo. 
At each location, monthly averages are formed from a 
minimum of 3 daily retrievals within each month: this 
threshold is easily exceeded by TOMS, but somewhat re- 
duces the number of AVHRR monthly averages that can 
be formed. At grid points lacking 3 daily retrievals in 
a particular month, the monthly average is constructed 
using linear interpolation between the monthly averages 
of the preceding and following months. Interpolation is 
applied only at grid points having at least 63 months 
with 3 daily retrievals or more (out of a maximum pos- 
sible 84 months); otherwise, the grid point is excluded 
from the analysis. Interpolation allows the correlation to 
include months when only one retrieval contains three or 
more daily images. (Nonetheless, our results are qual- 
itatively unchanged when these months are omitted.) 
The 63 month threshold is chosen to exclude from our 
analysis locations that are consistently devoid of data. 
We also correlate the retrievals limiting ourselves to the 
summer months (June to September, with 28 potential 
months over 7 years); the threshold is 21 months. Fi- 
nally, to compare the interannual variations of each re- 
trieval, monthly anomalies are constructed by forming a 
seven-year average of each of the 12 calendar months at 
every grid point and subtracting this seasonal cycle from 
the monthly averages. 

3. Comparison of Seasonal and 
Interannual Variations 

The spatial distribution of the TOMS AI and AVHRR 
AOT retrievals are shown in Fig. 2. The figure illustrates 
the various aerosol types contributing to the annual and 
summertime mean of each retrieval. Both show an opti- 
cally thick aerosol plume extending across the Atlantic 


from North Africa, which ground-based measurements 
identify as dust originating from the fringes of the Sahara 
during NH summer and additionally from the Sahel dur- 
ing winter [Car/son and Prospero , 1972; Chiapello et al . , 
1995; Li et a/., 1996; Chiapello et ai : 1999]. The plume 
extending westward off the coast of Southern Africa is 
attributed to carbonaceous aerosols formed by biomass 
burning during the dry season (e.g., [Hsu et al ., 1999]). 
The AVHRR retrieval shows an additional plume ex- 
tending eastward from North America into the North 
Atlantic, which is interpreted as sulfate aerosols created 
by industrial pollution [Husar et ai } 1997]. 

Based upon in situ observations (e.g. [Carlson and 
Prospero , 1972; Chiapello et a/., 1999], we define a ‘dust 
region,’ extending from 40. 5° W to 18.125°W longitude 
(near the African coast), and from 5.5°N to 25.5°N lat- 
itude, where soil dust is the dominant aerosol type. In 
this region, we expect that both the TOMS AI and 
the AVHRR AOT will record mainly the presence of 
soil dust, and that any potential agreement between 
the two retrievals will not be degraded by other aerosol 
types for which each retrieval has a different sensitivity. 
Figs. 3a and b show the seasonal cycle within the dust 
region. Longitudinal variations are small (not shown), 
except for a gradual decrease in amplitude to the west; 
in the figure, we average each retrieval longitudinally 
across the dust region. Near the equator, the aerosol 
load peaks during NH winter; as summer approaches and 
the ITCZ moves northward, the trade winds that trans- 
port the dust shift polew r ard as well, so that at north- 
ern latitudes, each retrieval peaks during the summer 
[Swap et a/., 1996]. A similar northward shift in dust 
transport during NH summer is indicated by Meteosat 
images [Moulin et a/., 1997a]. The shift is accompanied 
by seasonal changes in the land source region. Saharan 
dust is transported over the Cape Verde Islands (24° W, 
16° N) throughout the year, while the occurrence of east- 
erly Trade winds at the surface during the winter brings 
dust additionally from source regions in the Sahel [Chia- 
pello et a/., 1995]. The seasonal cycle is more pronounced 
in the TOMS AI retrieval, as indicated by Fig. 4, which 
shows the ratio of the AI and AOT as a function of cal- 
endar month, averaged along 15.5° N between the longi- 
tudes of Dakar and the Cape Verde Islands. Note also 
that the summertime maximum of the AI at northern 
latitudes is roughly twice the winter maximum that oc- 
curs to the south (Fig. 3a). In contrast, the maximum 
value of the AVHRR AOT shows little variation with 
latitude. 

These contrasts can be attributed to seasonal changes 
in the height of the aerosol layer. Aerosol transport off 
North Africa during the NH summer occurs in an ele- 
vated layer centered near 700 mb [Carlson and Prospero , 
1972], whereas wintertime transport is centered closer to 
the surface, near 900 mb [Chiapello et al 1995]. Con- 
version of the TOMS AI into an equivalent AOT at 
the 630 nm wavelength retrieved by AVHRR involves 
two steps. First, Fig. 1 of [Hsu et a/., 1999], which 
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is based upon calculations by [Torres et a/., 1998], is 
used together with the aerosol layer height to convert 
the AI into a 380 nm AOT. Second, this AOT is con- 
verted into a 630 nm value, using a theoretical model of 
AOT dependence upon wavelength and the Angstrom 
exponent parameter derived from sun photometer mea- 
surements. Based upon such measurements in regions 
where the aerosol load is dominated by dust, [Hsu et a/., 
1999] derive an Angstrom exponent that indicates the 
630 nm AOT is approximately the same as the 380 nm 
value. (More precisely, the latter is roughly 5% larger.) 
Thus, the most important factor in relating the AI to an 
AVHRR AOT is the aerosol height. Based upon Fig. 1 
of [Hsu et a/., 1999], the effect of raising the aerosol layer 
from its wintertime level near 900 mb to its summer el- 
evation of 700 mb would by itself increase the AI by a 
factor of just over two. This is consistent with the sea- 
sonal contrast of the ratio of the TOMS AI and AVHRR 
AOT shown in Fig. 4. This is also roughly the ratio of the 
summer and winter maxima of the AI (Fig. 3a). Thus, 
both retrievals suggest that the amount of dust trans- 
ported off Africa remains constant throughout the year, 
even as latitude of transport varies. Clearly, our attempt 
to infer seasonal variations in AOT from the TOMS AI 
by accounting for seasonal variations in the height of 
aerosol layer is limited by the uncertainty in this height; 
a more precise comparison would require detailed infor- 
mation on the vertical distribution of the aerosol trans- 
port. Nonetheless, the present comparison suggests that 
the seasonal cycles of AOT inferred from the two re- 
trievals are roughly consistent with each other; differ- 
ences in the two retrievals themselves can be attributed 
to seasonal variations in the height of the aerosol layer. 

Interannual variations in the two data sets, calculated 
by removing the seasonal cycle from the monthly aver- 
ages, have a spatial distribution resembling the time av- 
erage, as indicated by the standard deviation of the inter- 
annual variability, shown in Fig. 5. The standard devi- 
ation measures the typical amplitude of each retrieval s 
monthly anomalies; our correlation of these anomalies 
from the two data sets is meaningful to the extent that 
this amplitude (the ‘signal ) is larger than the uncer- 
tainty (the ‘noise 1 ) of each retrieval. The uncertainty 
of gridded daily images of the TOMS AI is estimated 
to be 0.2 [Hsu et a/., 1999], while that of the AVHRR 
AOT is estimated at less than 0.04 [Stowe et al ., 1997]. 
The uncertainty of the monthly averages and anoma- 
lies should be further reduced from the daily value by 
roughly the square-root of the number of independent 
observations available each month. Independent mea- 
surements of the dust aerosol amount are separated by 
roughly three days, according to Fig. 6a, which shows 
the lag autocorrelation of the TOMS AI falling to zero 
over this duration. We estimate the number of indepen- 
dent daily images available in each month by counting 
the number of successive images spaced three or more 
days apart. (For example, daily images available on the 
first, fourth, and seventh day of the month would be 
counted as three independent measurements, whereas 


images available on the first, second, and fifth would 
count only as two.) Figs. 6b and c show that within 
the dust region, roughly 4 AVHRR AOT images and 6 
TOMS AI images are independent each month. This is 
a conservative estimate of availability since images from 
adjacent days are not entirely redundant. In any case, 
the uncertainty of the monthly averages and anomalies 
can be estimated to be less than half of the uncertainty 
of daily images cited above. Over the NH subtropical 
Atlantic, this is exceeded by the standard deviation of 
each retrieval (Fig. 5). Thus, month to month varia- 
tions in the two data sets predominately record varia- 
tions in aerosol absorption and reflectance rather than 
measurement uncertainty, a necessary condition for our 
comparison of interannual variability. This dominance 
is especially large near the African coast, close to the 
source region, where the standard deviation is largest. 

Fig. 7a shows the month to month variations of each 
retrieval averaged over the dust region for the period 
of comparison. The two time series exhibit reasonable 
resemblance, as reflected by the correlation coefficient 
of 0.66. This agreement is partly due to agreement of 
the seasonal cycle in each data set (Fig. 7b); both re- 
trievals identify the same month corresponding to max- 
imum transport as a function of latitude (Fig. 3). Year 
to year variations can be compared by removing the sea- 
sonal cycle (Fig. 7b). As noted in the previous section, 
we assume that interannual variations in the height of 
the aerosol layer are of secondary importance to the AI, 
so that on interannual time scales, it can be compared 
directly to the AVHRR AOT. In the absence of the sea- 
sonal cycle, the correlation of the two retrievals is re- 
duced to 0.55 (Fig. 7c). The variance shared by the two 
data sets, equal to the square of the correlation coeffi- 
cient, is roughly 30 percent. 

When the two retrievals are compared at individual 
grid points, the removal of the seasonal cycle leads to a 
greater reduction in the correlation coefficient. Fig. 8a 
shows the spatial distribution of the correlation coeffi- 
cient prior to the removal of the seasonal cycle; locations 
without sufficient data to form at least 63 monthly av- 
erages are denoted by dots. The predominance of dust 
scattering over that by other aerosol types and the pre- 
dominance of the month to month variations of each re- 
trieval over measurement uncertainty, two assumptions 
upon which the validity of this comparison depends, 
are best satisfied within our dust region defined above. 
Nonetheless, the figure displays the correlation over a 
broader area than the dust region in order to encompass 
the substantial dust concentrations that are observed to 
the north [Moulin et ai y 1997a] and downwind toward 
Barbados [Li et a/., 1996] during certain months of the 
year. According to Fig. 8a, correlations that are high 
given the inclusion of the seasonal cycle, are generally 
reduced when the seasonal cycle is removed (Fig. 8b), 
Despite agreement on the seasonal cycle of dust loading, 
the TOMS AI and AVHRR AOT are not well correlated 
on interannual time scales. 
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The uncertainty in the TOMS AI retrieval of soil dust 
is potentially largest during the NH winter, when the 
aerosol layer is closest to the surface, resulting in a low 
value of the AI that could be adulterated by the presence 
of more weakly absorbing aerosols or clouds. In contrast , 
summertime transport of dust across the Atlantic occurs 
w f ell above the surface, where its A I signature is much 
stronger compared to that of other aerosol types. Fur- 
thermore, POLDER 1 estimates of Angstrom coefficients 
show evidence of biomass burning and soot aerosols in 
the southeast corner of the dust region during Novem- 
ber and December (at least for the year 1996; [Chiapello 
et a/., 2000]). To remove these effects upon our com- 
parison, which introduce uncertainty into the soil dust 
estimates, we repeat the correlation restricting ourselves 
to just the summer months (June to September), distin- 
guished by a darker line in Fig. 7. The correlation in 
the absence of the seasonal cycle, equal to 0.07, is ev- 
idently not improved compared to the value computed 
using all months, which is also the case at individual grid 
points, as shown in Fig. 8d. The disagreement between 
the monthly anomalies of the TOMS AI and AVHRR 
AOT is not due to the lower level of aerosol transport 
during NH winter, nor the presence of soot aerosols dur- 
ing this season. We also restricted the correlation to 
the period between Jan 1985 and Dec 1987, when the 
AVHRR retrievals w r ere provided almost exclusively by 
the NOAA-9 satellite; wdiile this minimized the effects 
of inter-satellite calibration, the correlations were not 
increased (results not shown). 

4. Estimation of Retrieval Monthly 
Averages 

The TOMS AI and AVHRR AOT provide consis- 
tent estimates of the seasonal cycle of the dust aerosol 
load over the NH subtropical Atlantic, if seasonal varia- 
tions in the height of transport are taken into account. 
However, the monthly anomalies share far less variance. 
Here, w r e suggest that the w^eak interannual correlation 
is at least partly the result of clouds that limit the avail- 
ability of daily images comprising the monthly averages 
of each retrieval (c.f. Fig. 1). Below', we show' that the 
two retrievals are highly correlated when the monthly 
averages and anomalies are constructed using only days 
when both retrievals are available. The large shared vari- 
ance suggests that both retrievals consistently and accu- 
rately measure the dust load of any particular scene. 
However, we also show that only the TOMS AI has a 
sufficient number of daily retrievals to estimate precisely 
the average for any particular month; the low interan- 
nual correlation of monthly anomalies comprised of all of 
the data (Fig. 8b) reflects in part the large uncertainty 
in the AVHRR AOT monthly anomalies. 

Fig. 9 show's the correlation of the two data sets, w hen 
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the monthly averages are restricted to days when both 
the TOMS AI and AVHRR AOT retrievals are available. 
Clouds are absent from these scenes, w f hich otherwise 
prevent retrieval by AVHRR and TOMS, the latter to a 
lesser extent. The monthly anomalies are correlated at 
a substantially higher level (Fig. 9b), compared to the 
anomalies formed using all the available days (Fig. 8b). 
The correlation is near 0.7 within the dust region, in- 
dicating that the tw'o retrievals share as much as half 
their variance. The correlation and shared variance are 
even higher when restricted to summer months (Fig. 9d). 
While the tw'o retrievals are highly correlated wdien their 
monthly averages are constructed from identical scenes, 
these monthly averages have a potentially large uncer- 
tainty, being constructed from only a small number of 
daily retrievals. Fig. 10 show's the number of common 
days available for each of the 84 months, averaged over 
the grid points within the dust region. Reflecting the 
relative paucity of AVHRR daily retrievals, this value is 
around 3 days with a standard deviation at individual 
grid points of just under 2. 

Next, w'e consider how many daily retrievals are 
needed to estimate precisely the monthly anomalies, and 
w'hether either retrieval satisfies this threshold. To es- 
timate the minimum number of daily retrievals needed 
to form an accurate monthly average, w^e correlate two 
versions of the TOMS AI monthly averages with the sea- 
sonal cycle removed from each. In one version, the full 
number of available days is used to form the monthly 
averages. In the second (or ‘degraded’) version, aver- 
ages are formed after removing observations at random. 
That is, the monthly averages in the degraded version 
are formed with less than the full amount of available 
retrievals. The question we ask is how many retrievals 
can be discarded before the two versions of the data set 
are no longer well correlated. Because the uncertainty 
in a degraded monthly average increases as days are dis- 
carded — along with the potential error — this calcu- 
lation provides an estimate of the minimum number of 
daily retrievals required to estimate accurately the aver- 
age of any particular month. 

Fig. 11 shows the results. For example, in the TO 
day 5 case (Fig. lie), an average for a particular month is 
formed using all the available daily retrievals at locations 
w here there are at least 10. In addition, a corresponding 
degraded version of the monthly average is formed by 
discarding daily retrievals at random so that each aver- 
age is based upon only 10 values. At locations w r here at 
least 63 months (out of 84) satisfy the 10 day thresh- 
old, remaining months w'ith less than 10 daily retrievals 
available are filled using linear interpolation between ad- 
jacent months. (Increasing the threshold from 63 to 74 
reduced the number of grid points available for compar- 
ison, but did not change the degree of correlation at the 
remaining locations.) To increase the number of months 
containing at least 10 days, the TOMS daily retrievals 
w'ere screened for clouds using a less stringent 25 percent 
reflectivity threshold. Within the dust region, this al- 


lows roughly 5 additional retrievals to contribute to any 
monthly average, so that 15 to 20 daily retrievals are 
available in an average month. This reflectivity thresh- 
old allows the two versions of the monthly anomalies 
to be correlated at a larger number of grid points, but 
again didn’t affect the agreement at locations contain- 
ing at least 63 months according to either threshold. At 
each grid point, the correlation was repeated 50 times, 
in each case with different days removed at random from 
the degraded version of the TOMS AI. The correlation in 
Fig. 11 is the average of the fifty calculations. (The aver- 
age correlation did not change very much as the number 
of calculations was increased past 10.) 

When the degraded version of the AI is limited to 
15 daily retrievals, it remains highly correlated with the 
monthly anomalies formed from the full TOMS AI data 
set, suggesting that 15 daily retrievals are sufficient to 
estimate accurately the monthly average. Similarly, the 
correlation is high when the degraded version is limited 
to only 10 days; apparently, the removal of 5 daily re- 
trievals causes only a small change in the monthly av- 
erage. However, as the number of daily retrievals is re- 
duced to fewer than 10 days, the two versions of the 
TOMS AI become increasingly decorrelated. If we in- 
sist that the original and degraded version of the TOMS 
AI monthly anomalies have correlation coefficients of at 
least 0.7 — so that at least half of the variance is shared 
— then at least 7 daily retrievals are needed to form reli- 
able estimates of the monthly average in the dust region. 

In deriving this threshold, we have not assumed or re- 
quired that the daily images comprising the monthly av- 
erages in each version are statistically independent (and 
thereby separated by three or more days, according to 
Fig. 6a). Thus, it is possible that the seven daily im- 
ages used to calculate the degraded monthly averages 
(Fig. lid) are not independent. In this case, fewer than 
seven daily images might be sufficient to derive accu- 
rate monthly averages if all the images are independent. 
Conversely, 7 retrievals bunched together early or late 
in the month will not be sufficient for a precise estimate 
of the monthly average, since some of the information 
contained within the 7 values is redundant. 

Fig. 1 shows that while the TOMS AI exceeds the 7 
day threshold throughout the dust region — even with 
the more stringent 15 percent cloud reflectivity crite- 
rion — the AVHRR AOT retrievals exceed this thresh- 
old only in the northern part of the region. This suggests 
that the poor correlation of the TOMS AI and AVHRR 
AOT in the absence of the seasonal cycle is at least partly 
due to the imprecise estimate of the AOT by AVHRR 
for any particular month, a result of the limited number 
of retrievals permitted by the presence of clouds. One 
means of increasing the number of retrievals available 
in any particular month is to reconstruct each data set 
using a coarser grid: that is, trade away spatial reso- 
lution in return for increased precision of the monthly 
averages. We found that we had to increase the size 
of the AVHRR grid box to 4° latitude by 5° longitude 


— resolution typical of a climate AGCM — in order to 
have at least 63 monthly averages of AVHRR compris- 
ing at least 10 days in the vicinity of the ITCZ. This is a 
small increase in retrieval availability given that the new 
grid was constructed by averaging together 20 AVHRR 
and 16 TOMS boxes from their original grids; presum- 
ably, locations of missing AVHRR data often occur adja- 
cent to each other, which obviates any advantage gained 
by averaging these points. Nonetheless, the correlation 
of the two retrievals after removal of the seasonal cy- 
cle is increased upon rebinning, as shown in Fig. 12b. 
This improvement could be interpreted to result from 
the greater availability of the AVHRR retrieval on the 
coarser grid. However, the correlation is not increased 
when restricted to the NH summer months (Fig. 12d). 
This is unexpected since the 7 day threshold can be de- 
rived by considering only summer months, and suggests 
that there are additional reasons for the disagreement of 
the two retrievals at interannual time scales. One pos- 
sibility is the presence of undetected sub-pixel clouds in 
the TOMS AI [Torres et a/., 1998], although we have not 
attempted to quantify this effect. 

5. Conclusions 

We compared two estimates of the seasonal cycle 
and interannual variability of soil dust aerosols, as in- 
ferred from the AVHRR AOT and TOMS AI satellite 
retrievals. Both data sets are attractive to aerosol mod- 
elers because of their multi-decadal span and global cov- 
erage. We undertook the comparison in order to iden- 
tify common features of variability that could be used 
to evaluate the dust distribution calculated by aerosol 
models, including AGCMs. Each of the retrievals is sen- 
sitive to a different aerosol property — reflectance in the 
case of AVHRR, and absorption for TOMS — and em- 
phasizes different aerosol types. We focused our atten- 
tion upon the NH subtropical Atlantic Ocean, where soil 
dust makes the largest contribution to the total aerosol 
load [Li et a/., 1996; Tegen et a/., 1997; Chiapello et a/., 
1999]. Within this region, we assume that soil dust dom- 
inates each retrieval, and that the variability associated 
with sulfate and carbonaceous aerosols is comparatively 
small. Both retrievals identify the same climatological 
month corresponding to maximum dust loading, includ- 
ing the variation of this month with latitude. Because 
the AI varies with the vertical distribution of the aerosol, 
the amplitude of its seasonal cycle cannot be directly 
compared with that of the AOT. However, when sea- 
sonal variations in the height of the aerosol layer are ac- 
counted for, the seasonal cycles of the dust load implied 
by each retrieval are comparable. 

When the seasonal cycle is removed from the monthly 
averages, the correlation between the two retrievals is 
sharply reduced. For this comparison, we followed the 
assumption of [Hsu et a/., 1999] that year to year varia- 
tions in the height of the aerosol layer are of secondary 
importance to the AI, which thus can be directly com- 
pared on interannual time scales to the AOT. In a sep- 
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arate study (to be reported upon elsewhere), we tested 
this assumption by comparing the AI and dust AOT 
produced by the Goddard Ozone Chemistry Aerosol Ra- 
diation and Transport model (Paul Ginoux, personal 
communication). This model calculates a realistic spa- 
tial distribution of dust aerosols by combining observed 
winds and precipitation with prescribed dust sources and 
sinks in a tracer transport model. We found that the 
monthly averages of the model A I and dust AOT are 
almost perfectly correlated, after the removal of the sea- 
sonal cycle. That is, while variations in aerosol height 
contribute strongly to the AI on daily and seasonal time 
scales, interannual variations of the AI are dominated by 
variations of the AOT. 

The poor interannual agreement is in contrast to the 
substantially larger correlation that results when the 
monthly anomalies are constructed using only days com- 
mon to the two retrievals. This consistency suggests 
that both the TOMS AI and AVHRR AOT accurately 
estimate the aerosol load of any particular scene. How- 
ever, because clouds limit the availability of the AVHRR 
AOT retrieval to a greater extent than the TOMS AI, its 
monthly averages have a potentially larger uncertainty. 
We suggest that this uncertainty is partly responsible 
for the poor interannual correlation, when the monthly 
anomalies are constructed using all available daily obser- 
vations. We estimated how many daily retrievals were 
needed to form a reliable estimate of the monthly aver- 
age by correlating two versions of the TOMS AI monthly 
anomalies. The anomalies in one version were con- 
structed using all the available data, typically 15 or more 
retrievals per month within the dust region. In the other 
version, a certain number of observations were discarded 
so that the monthly anomaly was estimated with greater 
uncertainty. When the degraded version was constructed 
with fewer than 7 to 10 retrievals per month, the corre- 
lation with the original TOMS AI was quite low. This 
suggests that estimates of interannual variations of soil 
dust transport over the Atlantic require roughly at least 
7 retrievals per month. An additional constraint results 
from the 3 day decorrelation time of dust, so that the 
7 retrievals must not represent consecutive days. Only 
the TOMS AI is able to exceed this threshold, due to its 
greater ability to retrieve the aerosol load in the presence 
of clouds. Whether it should be regarded as a more re- 
liable indicator of interannual variability depends upon 
the importance of undetected sub-pixel clouds [Torres 
et a/., 1998]. This is potentially a more serious limita- 
tion of TOMS, due to the instrument’s relatively large 
footprint, although we did not quantify this effect. In 
any case, both TOMS and AVHRR should be capable 
of accurately estimating the seasonal cycle, since seven 
years provide sufficient observations necessary to form 
the climatological mean of any calendar month. 

In order to increase the number of AVHRR AOT re- 
trievals available during any particular month, we re- 
binned the data onto a coarser grid. This resulted in an 
improved correlation with the TOMS AI, which we inter- 


preted as due to the increased precision of the AVHRR 
monthly averages. In addition, spatial averages of each 
retrieval over the dust region showed a higher interan- 
nual correlation, compared to the correlation at individ- 
ual grid points (c.f. Fig. 7c versus 8b). This suggests 
that there is intrinsic agreement between the TOMS 
AI and AVHRR AOT that is mainly obscured by the 
sampling uncertainty of the latter. That is, sub-pixel 
cloud contamination may be of lesser importance. How- 
ever, there was almost no agreement when the compar- 
ison of spatial averages was restricted to the summer 
months. Similarly, rebinning did not improve the sum- 
mertime correlation. Sub-pixel contamination by clouds 
may be particularly important during this season, when 
mesoscale and synoptic variations of cloud cover associ- 
ated with African waves is largest. 

The high correlation of scenes common to the TOMS 
AI and AVHRR AOT suggest that current estimates of 
interannual aerosol variability are limited not by lack of 
precision of either retrieval wdthin individual scenes, but 
rather by the insufficient number of AVHRR retrievals 
needed to estimate accurate monthly averages and thus 
month to month variability. However, the availability of 
the AVHRR AOT is not far below the seven-day thresh- 
old that we estimated; only a few more daily retrievals 
per month could make a large improvement. This might 
result if a more precise cloud screening algorithm could 
exclude fewer scenes. In addition, the uncertainty of the 
TOMS AI could be reduced by quantifying the effect of 
sub-pixel clouds. In the meantime, the disagreement be- 
tween the two retrievals at interannual time scales might 
be adjudicated using Meteosat retrievals [Moulin et a/., 
1997b], which despite their restricted areal coverage, in- 
clude our dust region. In addition, a more precise com- 
parison of the annual cycle of dust loading may be possi- 
ble using TOMS AOT retrievals that incorporate infor- 
mation on aerosol height more precisely and systemati- 
cally than we have done here [Torres et al. , 1998]. 

The limited interannual correlation may also indicate 
the presence in each data set of variability by other 
aerosol types. In addition to dust, AVHRR retrievals 
may include industrial sulfates, carbonaceous aerosols, 
and sea salt, while TOMS may include carbonaceous 
aerosols. We tried to minimize this difference by lim- 
iting our study to the NH subtropical Atlantic Ocean, 
where we expect soil dust to be the dominant source 
of aerosol variability [Chiapello et al ., 1999]. The in- 
troduction of polarization measurements in the retrieval 
will make available additional information, such as the 
particle index of refraction or Angstrom coefficient, that 
can identify the aerosol type [Mishchenko and Travis , 
1997; Chiapello et a/., 2000]. Similarly, two-channel 
AVHRR will replace assumptions about aerosol particle 
size with inferences derived from the additional channel 
[Mishchenko et a/., 2000]. 

Despite the limited variance shared by the two re- 
trievals at interannual time scales, the seasonal cycle of 
dust loading, for which both data sets agree, can be 
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used to constrain dust models. Certain global models 
of aerosols underestimate in comparison to the AVHRR 
AOT the summertime maximum of soil dust transported 
downwind from the Sahara [Tegen et a/., 1997; Perlwitz 
et al . , 2000]; for these models, the agreement of the sea- 
sonal cycle between the two retrievals is more imme- 
diately relevant than the uncertain estimates of inter- 
annual variability. Future retrieval products will lead 
to stronger constraints on model interannual behavior. 
Light Detection And Ranging (LIDAR) observations will 
provide complementary information on the vertical dis- 
tribution, which is crucial to the climate effect of aerosols 
[Miller and Tegen , 1999]. Together, these observations 
will lead to more confident predictions on the importance 
of aerosols in forcing climate and climate change. 
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Figure Captions 

Fig. 1: Average number of daily retrievals available per month from 1984 to 1990 for a) TOMS AI and b) A\ HRR 
AOT. 

Fig. 2: Climatological average of the TOMS AI and AVHRR AOT retrievals, a) Annual mean TOMS AI, b) annual 
mean AVHRR AOT, c) NH summer (June-September) mean TOMS AI, and d) NH summer (June-September) mean 

AVHRR AOT. 

Fig. 3: Seasonal cycle of the a) TOMS AI and b) AVHRR AOT, averaged zonally across the NH subtropical Atlantic 
from 40.5°W to 18.125°W longitude (near the African coast). 

Fig. 4: Seasonal cycle of the ratio of the TOMS AI and AVHRR AOT, averaged at 15.5°N between the longitudes 
of Dakar and Sal, Cape Verde Islands. 

Fig. 5: Standard deviation of the TOMS AI and AVHRR AOT monthly average retrievals with the seasonal cycle 
removed, a) entire year, TOMS AI, b) entire year AVHRR AOT, c) NH summer (June-September) TOMS AI, and 
d) NH summer (June-September) AVHRR AOT. 

Fig. 6: a) Autocorrelation as a function of lag for the TOMS AI. The 1988 AI daily retrievals were high-pass filtered 
to remove frequencies higher than 30 days. The autocorrelation is averaged over the dust region. (At any individual 
location, the function is quite similar.) b) Number of indepedent observations available each month for the TOMS 
AI, and c) AVHRR AOT. The number of indepedent observations was computed by counting successive daily images 
within each month separated by three or more days. 

Fig. 7: a) Spatial average within the dust region of the TOMS AI (solid) and AVHRR AOT (dotted) monthly 
averages. For ease of comparison, the spatial average of each data set has been normalized to have zero mean and 
unit standard deviation over the 84-month period, b) Seasonal cycle of the spatial average of each data set, c) 
Monthly anomalies of the spatial average computed by removing the seasonal cycle from the monthly averages in (a). 
Values corresponding to NH summer months (June-September) are denoted by thicker lines. 

Fig. 8: Correlation of TOMS AI and AVHRR AOT a) including the seasonal cycle, b) with the seasonal removed, 
c) restricted to the summer months, but including the seasonal cycle, d) restricted to the summer months, with the 
seasonal cycle removed. The 95%, 99%, and 99.9% confidence levels are a, b) 0.21, 0.28, 0.35, and c, d ) 0.37, 0.47, 
0.58. The correlation is computed at locations where at least 63 months have 3 or more daily retrievals; otherwise, 
the location is denoted by a dot. 

Fig. 9: As in Fig. 8, but with the monthly averages of each retrieval computed using only the days when both 
retrievals are available. 

Fig. 10: The number of days during any particular month when both retrievals are available at a grid point, averaged 
over the dust region. 

Fig. 11: Correlation (after removal of the seasonal cycle) of monthly averages of the TOMS AI with a degraded 
version of the TOMS AI where monthly averages are formed after discarding daily retrievals at random. Dady 
retrievals are discarded until the degraded version has a) 1, b) 3, c) 5, d) 7, e) 10, f) 15 retrievals available to form 
the monthly average. Correlations are computed at locations where at least 63 months have the minimum number 
of daily retrievals: e.g. in panel e , this minimum is 10. Otherwise, the location is marked in yellow. 

Fig. 12: As in Fig. 8, but with both the TOMS and AVHRR AOT rebinned onto a 4° latitude by 5° grid. Each 
monthly average is formed with at least 10 daily retrievals. 
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Figure 3. Seasonal cycle of the a) TOMS AI and b) AVHRR AOT, averaged zonally across the NH 
subtropical Atlantic from 40. 5° W to 18.125°W longitude (near the African coast). 



Figure 4. Seasonal cycle of t 
longitudes of Dakar and Sal, 1 


Ratio of TOMS (Al) to AVHRR (AOT) 



a) Annual Std. Dev. TOMS (84-90) 



0.0 0.1 0.2 0.3 0.4 0.50. 


Figure 5. Standard deviation of the TOMS AI and 
seasonal cycle removed, a) entire year, TOMS AI, b) 
September) TOMS AI, and d) NH summer (June-Septe 
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Figure 6. a) Autocorrelation as a function of lag for the TOMS AI. The 1988 AI daily retrievals were 
high-pass filtered to remove frequencies higher than 30 days. The autocorrelation is averaged over the dust 
region. (At anv individual location, the function is quite similar.) h) Number of indepedent observations 
available each month for the TOMS AI, and c) AVHR.R AOT. The number of indepedent observations was 
computed by counting successive daily images within each month separated by three oi more days. 
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o) Monthly timeseries of TOMS and AVHRR 
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b) Annual Cycle of TOMS and AVHRR 
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Figure 7. a) Spatial average within the dust region of the TOMS AI (solid) and AVHRR AOT (dotted) 
monthly averages. For ease of comparison, the spatial average of each data set has been normalized to have 
zero mean and unit standard deviation over the 84-month period, b) Seasonal cycle of the spatial average of 
each data set, c) Monthly anomalies of the spatial average computed by removing the seasonal cycle from 
the monthly averages in (a). Values corresponding to NH summer months (June-September) are denoted by 
thicker lines. 
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Figure 8. Correlation of TOMS AI and AVHRR AOT a) including the seasonal cycle, b) with the seasonal 
removed, c) restricted to the summer months, but including the seasonal cycle, d) restricted to the summer 
months, with the seasonal cycle removed. The 95%, 99%, and 99.9% confidence levels are a, b) 0.21, 0.28, 
0.35, and c, d ) 0.37, 0.47, 0.58. The correlation is computed at locations where at least 63 months have 3 
or more daily retrievals; otherwise, the location is denoted by a dot. 
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Figure 11. (Color figure) Correlation (after removal of the seasonal cycle) of monthly averages of the 
TOMS AI with a 'degraded 7 version of the TOMS AI where monthly averages are formed after discarding 
daily retrievals at random. Daily retrievals are discarded until the degraded version has a) 1, b) 3, c) 5, 
d) 7, e) 10, f) 15 retrievals available to form the monthly average. Correlations are computed at locations 
where at least 63 months have the minimum number of daily retrievals: e.g. in panel e , this minimum is 
10. Otherwise, the location is marked in yellow. 











22 




Correlation: 3 Days of TOMS 



Correlation: 5 Days of TOMS 




0.20 0.30 0.40 0.50 0.60 0,70 0.80 0.90 


0.65 




0.70 0.75 0.80 0.85 0.90 0.95 1.00 




Figure 11. (Color figure) Correlation (after removal of the seasonal cycle) of monthly averages of the 
TOMS AI with a ‘degraded’ version of the TOMS AI where monthly averages are formed after discarding 
daily retrievals at random. Daily retrievals are discarded until the degraded version has a) 1, b) 3, c) 5, 
d) 7, e) 10, f) 15 retrievals available to form the monthly average. Correlations are computed at locations 
where at least 63 months have the minimum number of daily retrievals: e.g. in panel e , this minimum is 
10. Otherwise, the location is marked in yellow. 
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Figure 12. As in Fig. 8, but with both the TOMS and AVHRR AOT rebinned onto a 4° latitude by 5 
grid. Each monthly average is formed with at least 10 daily retrievals. 



